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Abstract 

Bose systems, subject to the action of external random potentials, are considered. 
For describing the system properties, under the action of spatially random potentials 
of arbitrary strength, the stochastic mean-field approximation is employed. When 
the strength of disorder increases, the extended Bose-Einstein condensate fragments 
into spatially disconnected regions, forming a granular condensate. Increasing the 
strength of disorder even more transforms the granular condensate into the normal 
glass. The infiuence of time-dependent external potentials is also discussed. Fastly 
varying temporal potentials, to some extent, imitate the action of spatially random 
potentials. In particular, strong time-alternating potential can induce the appearance 
of a nonequilibrium granular condensate. 
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1 Introduction 



Physics of systems with a Bose-Einstein condensate in random media has been attracting 
attention for many years. At the beginning, the interest was concentrated on the behavior of 
hquid hehum in nanoporous media [1]. In recent years, the physics of dilute Bose gases has 
gained much interest [2-8] . Quasidisordered Bose gases have been reahzed experimentally by 
creating quasiperiodic optical lattices [9-11]. Disordered optical lattices have been realized 
by incorporating random impurities into an optical lattice [12]. There exists vast literature 
devoted to the theory of disordered optical lattices. A large number of references on the 
subject can be found in the recent review article [13]. 

A different type of disordered Bose systems corresponds to the systems without a periodic 
lattice potential, but which are subject to a spatially random potential. Experimentally, such 
random potentials are realized by means of optical speckles [14-17]. Physical properties of 
these Bose systems in an external random potential have been theoretically studied for the 
case of weak disorder [18-21] and for arbitrarily strong disorder [22,23]. 

The properties of uniform Bose systems and of those inside periodic lattice potentials 
are, of course, different [13] and may require differing theoretical methods of description. 
For example, the self-consistent mean-field approximation [24-29] describes well the uniform 
Bose systems with any strong atomic interactions. But the region of applicability of this 
approximation can be limited for Bose atoms in a lattice at zero temperature in the vicinity 
of the superfluid-insulator phase transition [30] . 

The properties of disordered Bose systems are also different, depending on whether ex- 
ternal random potentials have been imposed on an initially uniform or spatially periodic 
system [13]. One interesting feature, common for both types of systems, is that sufficiently 
strong disorder can lead to the occurrence of a new phase, called Bose glass. This state of 
matter, occurring in disordered lattices, was suggested by Fisher et al. [31] (see also the 
recent articles [32,33] and the review paper [13]). The Bose glass phase is often character- 
ized by the remaining presence of the condensate fraction (no > 0), but with the absence of 
superfluidity (ug = 0). 

The nature of the Bose glass phase is usually described as follows. The whole system 
fragments into the islands of the Bose-Einstein condensate localized in the deep wells of the 
random potential, while other regions, surrounding these islands, are filled by the normal 
fiuid, containing no condensate. This is why such a phase is also termed the granular 
condensate or the localized condensate. This phase has been observed in disordered optical 
lattices [34] and in nanoporous media filled by liquid helium [35]. Such a state of matter can 
also exist in Bose systems without lattices [13,36]. 

The present paper concerns the Bose systems without periodic lattices. In the absence of 
a random potential, such a system would be uniform, exhibiting at low temperatures Bose- 
Einstein condensation. The description of a Bosc-condcnscd system in external spatially 
random potentials of arbitrary strength can be accomplished by means of the stochastic 
mean-field approximation. We emphasize that describing such random systems requires 
a special caution, since perturbation theory is not always applicable to them. This, e.g., 
concerns the weak-disorder expansion, which can become invalid for random Bose-condensed 
systems. For these systems, the perturbation theory with respect to the coupling parameter 
can also fail. Therefore the use of more refined approaches, such as the stochastic mean-field 
approximation, is of great importance for the correct description of random systems. 
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After presenting the main characteristics of Bose-condensed systems in spatial random 
potentials, we turn to the problem of these systems in temporal alternating potentials. We 
show that the action of temporal external fields can induce in a Bose system the consequences 
analogous to those produced by spatially random potentials. In particular, a nonequilihrium 
granular condensate can arise under sufficiently strong external fields. Very strong alternat- 
ing fields will destroy the granular condensate, transforming the whole system into a normal 
turbulent fluid. 

Throughout the paper, the natural system of units is employed, where %— 1 and — ^■ 



2 Bose Systems in Random Potentials 



We consider a dilute Bose system composed of atoms interacting through the local potential 

a. 



$(r) = ^q5{v) , $0 = 47r 



m 



(1) 



in which is the scattering length and m, atomic mass. The system is subject to a spatially 
random external potential ^(r). So that the energy Hamiltonian is 



H 



^t(r) 



2m 



+ e(r) 



■^'''(r)'^'''(r)'^(r)'^(r) dr , 



(2) 



where is the Bose field operator. 

Without the loss of generality, the random potential can be taken as zero-centered, such 
that its stochastic averaging gives 



< iiy) > = . 

The stochastic correlation function 

R(r - r') = « e(r)e(r') » 
is assumed to be real and symmetric, 

R*{y) = R{-r) = R{r) . 



(3) 



(4) 



(5) 



The presence of a Bose-Einstein condensate implies the spontaneous gauge symmetry 
breaking [8] , which is standardly realized through the Bogolubov shift [37] of the field oper- 
ator 

iPir) ^ V^(r) = 77(r)+Vi(r). (6) 
Here ?7(r) is the condensate wave function normalized to the number of condensed atoms. 

No ^ J \vir)\' dr , (7) 

and "01 (r) is the Bose field operator of uncondensed atoms defining the number operator 

TVi = y ^I(r)^i(r) dr (8) 
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of uncondensed atoms. The field variables ri{r) and ^i(r) are treated as two independent 
variables, orthogonal to each other, 

j ri*{r)i;i{r) dr = . (9) 

Under the Bogolubov shift (6), in order to make the theory self-consistent, it is necessary, 
as is proved in Refs. [24-29], to introduce the grand Hamiltonian 

H^H-lioNo-iiiN, , (10) 

in which H is the energy Hamiltonian (2), with the Bogolubov-shifted field operators (6), 
and /Iq and /ii are the Lagrange multipliers guaranteeing the theory self-consistency. 

In the presence of random fields, there exist two types of averaging. One type is the 
stochastic averaging, denoted as <^ ... 3>, which characterizes the averaging over the dis- 
tribution of random potentials. And there is, as usual, the quantum statistical averaging, 
which for an operator A is defined as 

<i>^= Trpi, (11) 
where p is a statistical operator. The latter, for an equilibrium system, is 

with (5 = 1/T being inverse temperature. The total average of an operator A, 

<i> = <Trpi> (13) 

includes both, the quantum and stochastic averaging procedures. The grand thermodynamic 
potential, corresponding to the frozen disorder, is 

Q = -T < In Tr e'^" > . (14) 

One should also keep in mind the quantum-number conservation condition 

<V'i(r) >= 0. (15) 

To satisfy the latter, the grand Hamiltonian (10) should be complimented by one more term 
guaranteeing the absence in H of the terms hnear in ■^i (r) , as is shown in Ref . [29] . Here we 

do not add explicitly such a linear killer, since for a uniform system or for a system uniform 
on average, linear in V'i(r) terms do not arise and condition (15) is automatically satisfied 
[24-29]. 

For the zero-centered random potential, for which property (3) holds, the system can be 
treated as uniform on average. Then one can set 

^(r) = Vp^ (^po = ^) , (16) 

where V is the system volume. The field operator of uncondensed atoms can be expanded 
over plane waves as 
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The field operators in the momentum representation, a^, define the momentum distribution 

Uk = < alak > (18) 

and the anomalous average 

CTk = < CLka-k > ■ (19) 
The density of uncondensed atoms is 

Pi = ^E^^' (20) 



and the total anomalous average is 



^i = ^E^^- (21) 

^ k^O 



With the total number of particles N, the system average density is 

p = - = po + Pi • (22) 
The density of uncondensed atoms (20) consists of two terms, 

Pi = Piv + Pg , (23) 

among which p^r is due to thermal fluctuations and interactions, while pq is caused by the 
random potential and defines the density of a glassy component or, briefiy speaking, the 
glassy density 

pG=y / < I < V'i(r) >H r > dv . (24) 
The general definition of the superfiuid density [3,29] results in the expression 



< P2 > 

in which P is the momentum operator 

P = j iP{t) {-iV) i){T) dr = Jipi{r) (-iV) t/'i(r) dr . 

It is convenient to define the dimensionless atomic fractions. Thus, the condensate frac- 
tion 

no = — ^l-ni (26) 
P 

is expressed through the uncondensed- atom fraction 

^1 = = ^ E ^fc ■ (27) 
P k^o 
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Equation (21) gives the anomalous fraction 

--7 = ^E-- (28) 

The uncondensed-atom fraction (27) is the sum 

ni = UN + no (29) 
of the normal fraction un = Pn/ P and the glassy fraction 

nG = y = ;^ I ^\< ^i(r) r » . (30) 
Finally, from Eq. (25) we have the superfluid fraction 

= 1 - 3^ > (31) 

in which 

< P2 > 

Q = (32) 

is the dissipated heat per atom. With expansion (17), the glassy fraction can be written as 

= E < > ' (33) 

k 

where 

ak =< ttk >H ■ (34) 

Together with expanding over plane waves the field operator of uncondensed atoms, as 
in Eq. (17), let us expand the random potential as 

e(r) = ^ E CkC^-' , ^' = iv I ■ ^^^^ 

The Fourier transformation of the correlation function (4) is 

^(r) ' " / ■ (36) 

Then, Fourier-transforming Eq. (4), we get 

<^ Ck^P > = hpRk . (37) 

Passing to the momentum representation, we accomplish in Hamiltonian (2) the Bogol- 
ubov shift (6) and substitute there the Fourier transforms (17) and (35). As a result, the 
grand Hamiltonian (10) acquires the form 

4 

H=J2H^''^ + He,t, (38) 

n=0 



6 



in which 



^^^^ ^^T. E ^oaWpak-gap+g , (39) 

9 fe,p(^0) 



with the part 

Hext = PoCoVv + y/p^ Yl + Ck^k) + i E , (40) 

which is due to the action of the external random potential. 



3 Stochastic Mean-Field Approximation 

To treat Hamiltonian (38), an approximation is needed. The third- and fourth-order terms 
in the operators ak, occurring in Eq. (39), can be simplified by means of the Hartree-Fock- 
Bogolubov approximation, as in Refs. [24-29]. However, the interaction of the random 
potential with atoms, described by part (40), cannot be treated in the simple mean-field 
procedure, since 

<ajt> = <Cjfc> = 0, (41) 

which would kill the last term in Eq. (40). Such a case would correspond to considering 
asymptotically weak disorder. To treat the last term in (40), a more delicate decoupling 
procedure is required. For this purpose, we shall employ the stochastic mean-field approxi- 
mation suggested and used earlier for other physical systems [38-43] . Following the idea of 
this approximation, we simplify the last term of Eq. (40) writing 

alap^k-p = [al^^p + o^lap - alapj ^k-p , (42) 
where is the quantum average (34). Under equation (42), we have 

< alup^k-p > = < alap^k-p > • (43) 

Approximations (42) and (43) allow for the consideration of arbitrarily strong disorder 
strengths [22,23]. 

The next nontrivial procedure is the diagonalization of Hamiltonian (38) by means of the 
nonuniform nonlinear canonical transformation 

ak = Ukbk + fe&Lfc + Wk(pk , (44) 

in which Uk, Vk, Wk, and ipk are to be defined so that the resulting Hamiltonian be diagonal 
in the operators bk, such that 

<bk>H = < bkbp >H-^ . (45) 
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Then the transformation (44) gives 

oik = < Ofe >i? = Wk(pk ■ 

Accomphshing the diagonahzation, we find 



u 



2 _ + £fc 



mc 

2ir 



1 



where 



2m 



2 2 I / 2\^ 
^fe = £fc + (^c j 



and £fe is the Bogolubov spectrum 



with the sound velocity given by the equation 

mc^ = (no + cr)p$o 
The random variable (pk satisfies the integral equation 

Vk = VPo Cfc 



1 ^k—p^p 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



of the Fredholm type. The nonuniform nonlinear transformation (44), with Eqs. (47) to 
(51), results in the grand Hamiltonian 



in which 



k 



Eb^ \ Y.i^k - ^k) - [l - rio(l + (7) + ^ (l - + a^) 



(52) 



With the diagonal Hamiltonian (52), it is straightforward to find the momentum distri- 
bution (18), 

(53) 



^k 

Uk — coth 



e 

2r 



and the anomalous average 



mc 



coth + < lafcl' > 



The random variable (46), in view of Eqs. (47), is 

Vk 



Oik = 



8 



(54) 



(55) 



^Prom here 



,2 < Ifkl^ > 



The fraction of uncondensed atoms (27) is represented as sum (29), in which 

dk 

coth ( ) - 1 - 



2p 

while the glassy fraction (33) becomes 



— f — coth ( 
2p J Lfffe \2T) 



(27r)3 ' 



(57) 



_ 1 /■ <Mi> /.ox 
" p y (a;fe + mc2)2 (2^ " ^^^^ 

The anomalous fraction (28) is also the sum 

a ^ aN + nG , (59) 

where 

and the no is the same as in Eq. (58). The terms un and are caused by thermal 
fluctuations and interactions, while the glassy fraction hq is due to the action of the random 
field. 

In the superfluid fraction (31), the dissipated heat (32) reads as 

Q^Qn + Qg: (61) 

with 

1 /■ /c^ dk 

^""^8^ J smh'{ek/2T) (2^ ^ ' 

being due to interactions and finite temperature, while 

/^^T^^cothfii) * (63) 



2mp J Ekicok + mc^) \2Tj (27r)3 

is the heat dissipated by the glassy fraction. 

Note that deriving the above formulas we have used the form of the Lagrange multiplier 

Po = {l + ni + a)p$o , (64) 

obtained by minimizing the grand potential (14) with respect to the number of condensed 
atoms A^O) a-nd the multiplier 

/ii = (1 + ni - a)p^o , (65) 

found from the condition of the condensate existence [13], which is equivalent to the require- 
ment that the spectrum of collective excitations be gapless. 
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4 Failure of Weak-Disorder Perturbation Theory 



Uniform Bose-condensed systems under the actfon of random potentials are usually con- 
sidered in the case of asymptotically weak interactions and weak disorder. Exceptions are 
the articles [22,23], where the stochastic mean-field approximation [38-43] was employed 
allowing for the description of systems with arbitrarily strong atomic interactions and arbi- 
trarily strong disorder. It is necessary to stress that the perturbation theory with respect 
to the disorder strength may be inapplicable to the Bose systems in random potentials. In 
the present section, we show that the weak-disorder perturbation expansion may lead to 
incorrect results. 

Let us consider the last term 

Hran = \[Nq (66) 

in Hamiltonian (52), explicitly describing the influence of the random potential. The related 
contribution to the internal energy is 

Eran = ^ H^dn ^ ) (67) 

which results in 

Eran = < > ^ = < V^O > • (68) 

Suppose, one considers weak disorder and assumes that the weak-disorder perturbation the- 
ory should be valid. Then, the first term in Eq. (51) can be treated as the zero-order 
approximation 

Vl^^-y/p'o^k. (69) 
Iterating with Eq. (69) the second term of Eq. (51) gives the first-order approximation 

(1) / — t Cfc-pCp /-r>N 
-VPoCk TT- 2^ 2 ■ (^0) 

V p u!p + mc' 

Employing approximation (70) makes it easy to derive all results obtained by other re- 
searchers using the weak-disorder perturbation theory (see Ref. [22]). For example, the 
internal energy (68), due to random fields, becomes 

(1) ^_Noy. « m » 

Using here the correlation formula (37) and passing to integration yields 

£;(!) ^ _ / ^o-^P ^ (721 

J {up + mc2) (27r)3 ■ ^ ' 

It is exactly this expression (72) that has been rederived by all authors who have used the 
weak-disorder perturbation theory. Equation (72) tells us that random fields diminish the 
internal energy. 

But the result is very different, if no perturbation theory has been involved. According 
to Eqs. (34) and (41), we have 

<afe> = <afc> = 0. (73) 
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Then, from Eq. (46), it follows that 



<(^jt>= 0. (74) 
Consequently, the random-field contribution to the internal energy (68) is exactly zero, 

Eran = < Hran > = . (75) 

This conclusion sounds reasonable, if one remembers that the zero-centered potential, with 
the zero mean (3), is considered. Thus, the perturbation-theory formula (72) is in contradic- 
tion with the exact expression (75), hence, Eq. (72) docs not seem to be correct. The same 
concerns other formulas that have been derived by means of the weak-disorder expansion. 
Such formulas do not seem to be reliable. The reason for the weak-disorder perturbation 
theory failure is the noncommutativity, in some cases, of the operations of expanding over 
the random fields and of averaging over these fields. 



5 Random Potential with Local Correlations 

The consideration of the previous sections is general, beeing applicable to any spatially 
random potentials. To specify the problem, it is necessary to concretize the type of the 
random potential. When the correlation length of the correlation function (4) is sufficiently 
short, at least much shorter than the healing length, then the random potential can be 
modelled by the Gaussian white noise, with the local correlation function 

Rir) = i?o5(r) . (76) 

Then the Fourier transform (36) gives Rk = Rq and Eq. (37) becomes 

< CkCp > = SkpRo . (77) 

For numerical investigation, it is convenient to introduce dimensionless variables. The 
interaction strength is characterized by the gas parameter 

7 = p"^as . (78) 



The dimensionless temperature is defined as 



For the dimensionless sound velocity, we have 

mc 



pl/3 ■ 



(80) 



Disorder is known to possess the property of localizing atomic motion inside the regions 
of the characteristic localization length, which can be estimated as the Larkin length [44] 

_ 47r 

^''"^ = 1^ ■ ^^^^ 



11 



The strength of disorder can be described by the disorder parameter 

C ^ f {p'"a = l) , (82) 

Hoc 

where a is the mean interatomic distance. 

If the random potential ^(r) is hmited by a finite amplitude Vr, such that 

\m<VR 

for all r pertaining to the considered system, and if the correlation function (4) has a finite 
correlation length Iji, then the parameter Rq in Eqs. (76) and (77) can be represented as 

Ro ^Vlll. (83) 

In that case, the locahzation length (81) becomes 

47r 

which is close to the expression given in Ref. [17]. Returning to the local correlation function 
(76) implies that /i? — > and Vr — > cxd, so that product (83) be finite. 

When disorder is asymptotically weak, so that i?o — > 0, then Zjoc extends to the size of the 
whole system. The latter then represents the standard Bose-condensed system with extended 
condensate. With growing disorder, when the localization length becomes much smaller than 
the system linear size L, but when lioc is yet much larger than the mean interatomic distance 
a, that is, when 

a < hoc < L , (85) 

then the Bose-Einstein condensate fragments into multiple pieces separated by the normal 
phase with no gauge symmetry breaking [34-36,45-47]. This type of matter is termed the 
Bose glass or the granular condensate. Finally, when disorder is so strong that koc ~ a, then 
no condensate is possible, since atoms are localized separately, each of them being trapped 
in a deep randomly located well of the random potential. The latter phase forms the normal 
glassy matter. Such a phase was observed [14] in a strong random potential created by laser 
speckles. At zero temperature, the phase portrait, in the variables of the gas parameter 
(78) and disorder parameter (82), should look as in Fig. 1. The granular condensate starts 
appearing when li^c ~ (10 — 100)a, hence C, ~ 0.01 — 0.1. And the granular condensate 
transforms into the normal glass when li^c ~ a, so that C ~ 1- The phase transition between 
the extended condensate and granular condensate is continuous, while that between the 
granular condensate and the normal glass is of first order. 

The granular condensate is a phase that cannot already be treated as uniform on av- 
erage, as it is done for the extended condensate. The granular condensate is a principally 
nonuniform condensate, which requires a separate consideration taking into account the spa- 
tial nonuniformity. Assuming that the system is uniform on average does not distinguish 
between the extended and granular condensates, but this description covers all the region, 
where any condensate is possible. However, the phase transition between the system with 
a condensate and the normal glass can be described. Keeping this in mind, we follow the 
consideration of the previous sections. 
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(84) 



First, we need to solve the integral equation (51). A good approximate solution to this 
equation is 



<Pk = r^— f ■ (86) 

This random variable enters the above equations, such as Eqs. (53), (54), (55), (58), and 
(63), in the form ^ Iv'fcP To find the latter, we employ the self-similar approximation 
theory [48-50] in the variant involving the self-similar factor approximants [51-54]. As a 
result, we obtain 

2 _ UqRqS^!'^ 



Then the glassy fraction (58) is 



= 7,4/7(, _ ■ (88) 



Note that the disorder parameter (82) naturally appears in Eqs. (87) and (88). Therefore 
these equations can be considered as defining the disorder parameter (82) as such. 

In dimensionless variables (78) to (80), the dimensional equation (50) for the sound 
velocity becomes 

s2 ^ 4^^(i _ni + a) . (89) 
Due to Eqs. (29) and (59), this can be rewritten as 

- 47r7(l - + (7;v) ■ (90) 
For the normal fraction (57), we have 



coth 



2 ^ 

S X 



dx\ . (91) 



2t 

According to expressions (22) and (23), the atomic fractions are normahzed as 

rio + ?^jv + = 1 • (92) 
Using this, the glassy fraction (88) can be represented as 

(l-njv)C /n^x 

The superfiuid fraction (31) contains the dissipated heat (61). The part (63) of the 
dissipated heat, due to the dissipation on the glassy fraction, in the case of the white noise, 
contains Eq. (87). Then, integral (63) diverges, but can be regularized [22]. So that, finally, 
for the superfiuid fraction, we find 

3/2 



4 c5 roo X { — 1 ) dx 

n^^l - - no 7= / ; . (94) 

The anomalous fraction (60), because of the local interaction (1), also diverges and 
requires to be regularized [13,24-29,55]. In this regularization, the asymptotic properties of 
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function (60), with respect to temperature, should remain correctly defined. Thus, in the 
limit of zero temperature, using the dimensional regularization, we have [13,24,27-29] 



(7jv ~ — ^y^T-fno 



(95) 



While, when temperature tends to the critical point 



tc = 3.312498 , 



(96) 



then the correct limit of Eq. (60) is 



St 



{t ^ tc) . 



(97) 




The critical point (96) is the same as for the ideal Bose gas, as it has to be for a mean-field 
approximation [29]. 

It is important to stress that taking account of the anomalous averages is principally 
necessary. It would be absolutely incorrect to use the Shohno trick [56] by omitting the 
anomalous averages, as one often does. This is mathematically wrong, since the anomalous 
averages at low temperature can be larger than the normal fraction rijv, and they are of 
the same order at finite temperatures below the critical point [23,57]. Keeping what is 
of the same order, but neglecting what can even be larger, cannot be called a reasonable 
approximation. It is also straightforward to show [24,29] that omitting the anomalous average 
renders the consideration not self- consistent and the system unstable. The origin of the 
resulting inconsistency is very easy to understand [29]. The existence of the anomalous 
average is due to the gauge symmetry breaking. The latter is also the cause of the Bose- 
Einstein condensate existence. Hence both, the anomalous average and the condensate, 
either exist together or do not arise at all. If one wishes to omit the anomalous average, 
then, to be self-consistent, one must neglect the condensate existence. Or, when the latter 
is assumed, one has to retain the anomalous average as well. The omission of the anomalous 
average, in addition to breaking the system stability and making the thermodynamics not 
self-consistent, also distorts the phase transition order, provoking a first-order transition. 
The latter is, of course, incorrect, since the Bose-Einstein condensation is the second-order 
phase transiiton, irrespectively to the interaction strength [29]. 

The correct asymptotic behavior of the anomalous average in the vicinity of the critical 
point tc, as in Eq. (97), guarantees the second-order phase transition for any value of the gas 
parameter [13,28,29]. This can be demonstrated by direct numerical calculations [13,28] and 
also by expanding the quantities of interest in powers of the relative temperature deviation 



(98) 



in the vicinity of tc- Then we obtain the dimensionless sound velocity 




(99) 



the condensate fraction 



3 



8 ' 



(100) 
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the anomalous average 



and the superfluid fraction 



(Tat 



3 3 

T + - 

2 8 



1 + 



3 

' 2 



1.741 



(101) 



(102) 



^From these asymptotic expressions, the second-order phase transition is evident. 

Trying to interpolate the anomalous average between the asymptotic hmits (95) and (97), 
we can reorganize Eq. (60) to the identical form 



On 



= -1 / 

2p J 



mc 



dk 



Sk (27r) 



1 

2p 



mc 



coth ( ^ ) - 1 



(103) 



Employing for the first term of Eq. (103) the dimensional regularization and passing to the 
dimensionless quantities (78) to (80), we get [24,27-29] 



2^2 



On — 



(2 



2 s"" r 
W Jo 



(x/r+^-i) 



1/2 



coth 



- 1 



dx . 



(104) 



This expression gives the correct low-temperature limit (95). However, the critical behavior 
(97) becomes disturbed. This disturbance, as numerical calculations [28] show, is not essen- 
tial for the gas parameter in the range < 7 < 0.3, for which the phase transition remains 
continuous. For larger 7, the use of expression (104) would result in the discontinuity at the 
critical point tc- Therefore, for 7 > 0.3, form (104) can be used below the critical region, 
while close to Tc, it has to be replaced by the correct limit (97) guaranteeing the continuous 
phase transition [28]. 

At zero temperature, the normal fraction (91) and the superfluid fraction (94) reduce to 



riN = 



371 



2 ' 



no . 



(105) 



These equations, together with Eqs. (90), (92), (93), and (95), have been analysed in detail 
[22,23], demonstrating a first-order phase transition on the line C, = C{l) between the system 
with a condensate and the normal glass with no condensate. 



6 Simple Model with Quenched Disorder 

In Sec. 4, we showed that the weak-disorder perturbation theory fails in the description of the 
Bose-condensed system in spatially random potentials. Therefore, for correctly describing 
such systems, more refined approaches are necessary, for instance, as the stochastic mean- 
field approximation of Sec. 3, which allows for the consideration of arbitrarily strong disorder. 

Here we show that, in the presence of disorder, perturbation theory with respect to atomic 
interactions can also fail. We demonstrate this by a simple model with quenched disorder, 
which allows for an explicit illustration of the point where the weak-coupling perturbation 
theory fails. 
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Let us consider the model Hamiltonian 

H^{l + i)^'' + g^^ (106) 

of the so-called zero-dimensional </7^-theory in the presence of quenched disorder. Here the 
variable </? e (—00, 00) imitates the Bose condensate function. The coupling parameter 

(7 > characterizes the interaction strength. And the random variable describes an external 
random potential, with a distribution p(4) . As is usual, we assume that the random potential 
is zero-centered, such that 

/oo 
= . (107) 
-00 

Its dispersion is given by 

/oo 
epiO ■ (108) 
-00 

The often used examples of the random-variable distribution are the Gaussian distribu- 
tion 

^ (- S) ^'''^ 

and the uniform distribution 

PuiO = ^ 0{D - \^\) . (110) 
Their dispersions, respectively, are 

The partition function of the system in the frozen random field is 

Z{g,i)^— / e-^'d^. (Ill) 



'■n 

The free energy of the system with quenched disorder writes as 

/oo 
p(01nZ(^,0c?e- (112) 
-00 

In the limit of no disorder, when — > 0, then, in both cases (109) and (110), one has 
p(0 ~^ <^(0- ^ result, 

/(^)^-lnZ(^,0) (A^^O) 

with 

Z{g, 0) = ^ / exp (-(^^ - g^p'^) d(p . 

The case of no disorder allows for the use of the weak-coupling perturbation theory, provided 
the latter is complimented by a resummation procedure [58,59]. 
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However, for finite disorder, the situation is more complicated. The partition function 
(111), with Hamiltonian (106), reads as 

Z{g,^)^-= / exp{-(l + 0¥^'-^<^'} (113) 
Expanding here the integrand in powers of g and integrating gives the series 

oo 

^(^,0 = E^n(05" (114) 
n=0 

with the coefficients 

(-l)"r(2n + l/2) 

v^r(n + i)(i + 02"+^/' ■ ^ ' 

Substituting series (114) into the free energy (112) shows that the latter diverges for all 
n = 0,1,2,..., in the case of the Gaussian distribution , because of the pole at ^ = — 1 in 
coefficient (115). In the case of the uniform distribution (110), integral (112) is finite only for 
weak disorder, for which D < 1 and < 1/3. But for any stronger disorder, with -D > 1, 
the free energy is not defined, since Eq. (112) diverges. This shows that perturbation theory 
with respect to the coupling parameter, generally speaking, fails for the disordered system. 

In order to develop a perturbation theory, some special tricks arc necessary. For example, 
the disorder strength and the interaction strength could be treated not as independent 
quantities, but as being related through the ratio 

A.|, (116) 

which is assumed to be fixed [60]. Then the free energy (112), for the Gaussian distribution 
(109), takes the form 

This expression can be expanded in powers of g, yielding in the fc-order 

/,(5,A) = X:an(%". (118) 

n=l 

The coefficients a„(A) are not easy to define, and for < A < 1 they were found to be 
represented as the sum 

an(A) = 6n(A)+c„(A), (119) 

whose terms for n ^ 1 are [60] 

bM = (-1)"+^ ^ ^ (1 - A)"-V2 , 

V2 TT 

-jj^ exp (^—^^/n + aj cos( (^/J^y/n + (3^ , 
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where 

2.959237 3.184867 In 2 ^ 37r 

Although the found coefficients a„ = a„(A) are vahd only for n ^ 1, formally using them 
for all n>l, and fixing A = 0.1, gives 

ai = 0.854049 , 02 = -3.07481 , 03 = 22.1387 , 

04 = -239.098, 05 = 3443, 06 = -61974.1, 

and so on. These coefficients grow very fast, making series (118) strongly divergent. The 
series were shown [60] to be so widely divergent that they could not be Borel summed. 

But the divergent series (118) can be summed employing the self-similar approximation 
theory [48-50]. We have accomplished the resummation procedure by using the self-similar 
factor approximants [51-54]. The latter have the form 

r,{g,X)^a^{X)gl[{l + Agr, (120) 
1=1 



where 



k/2, /c = 2,4, ... 

{k + l)/2, A; = 3,5,... 



The coefficients and rii are defined by re-expanding form (120) and equating it to series 
(118), which yields 

J^riiA^^Br, (n = l,2,...,fc), 

1=1 

(n - 1)! 9^0 dg'"- [ ai{X)g 

and Ai = 1 for the odd Nk = (A; + l)/2. For the fixed A = 0.1 and different g, we calculated 
the factor approximants (120) up to the order k = 10. For g = 0.1, we get f^^ = 0.0662, 
with the relative error 22% as compared to the exact value 0.0542. For = 1, we have 
/iQ = 0.229, with an error 19%, as compared to the exact value 0.251. And for g — 10, we 
find fiQ = 0.712, deviating from the exact value 0.642 by 11%. These approximants should 
be considered as quite good, if one remembers that expressions (119) have been defined 
only for n ^ 1, but used for all n = 0, 1, 2, ... , 10. Therefore the resulting errors rather 
characterize the inaccuracy of these coefficients (119), but not merely the accuracy of the 
factor approximants (120). 

The main message of this Section is that, under disorder, the perturbation theory with 
respect to atomic interactions can become inapplicable. For the Gaussian distribution of the 
disorder potential, this perturbation theory fails for any strength of disorder, and for the 
uniform disorder distribution, the perturbation theory becomes invalid at sufficiently strong 
disorder. To realize the weak-coupling perturbation theory, it is necessary to invoke some 
tricks, like fixing a relation between the disorder strength and the interaction strength. How- 
ever, resorting to such tricks looks too artificial, since in reality the disorder and interaction 
strength are independent characteristics. 
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7 Creation of Nonequilibrium Granular Condensate 



In the previous Sections, we have considered equihbrium systems with Bose-Einstein con- 
densate, which, by applying spatially random fields, could be transformed into a Bose-glass 
type phase with a granular condensate. In the present Section, we show that a Bose system, 
subject to the action of an external time- alternating potential enjoys the properties to some 
extent analogous to those of a system in a spatially random potential. In particular, in such 
a nonequilibrium system a nonequilibrium granular condensate can be created. 
Let us consider a Hamiltonian 

H{t)^H + V{t) , (121) 

in which the first term 



f ^t(r) - ^ + U{r) ^{r) dr + ^ f 7/'t(r)^t(r)^(r)V'(r) dr (122) 

describes Bose atoms in a trapping potential U (r) , and the second term 

V{t) = [ ilj\r)V{r,t)ij{r) dr (123) 



corresponds to an external alternating potential. Again, for brevity we do not write explicitly 
the time dependence of the field operators ip{r) = ip{r,t). 

The external potential V{r,t) is assumed to vary in time so that its characteristic varia- 
tion time tyar is much larger than the local-equilibrium time tioc, but much shorter than the 
observation time, 

tloc tvar tabs • (-^^'^) 

As an estimate of the local-equilibrium time, we can take the values of the parameters typical 
of the dilute trapped gases ^^Rb and ^'^Na, say, m ~ 10~^^ g, P ~ 10^^ cm~^, and ~ 10~^ 
cm. The local-equilibrium time is defined [3] as 

in 

tioc = ■ (125) 

pas 

With the above parameters, this gives tioc ~ 10~^ s. 

Under condition (124), it is possible to define the local-equilibrium free energy 

F{t) = -T{t) In Tre-^(*)-^(*) , (126) 

in which T{t) = 1/ j3{t) is the local-equilibrium temperature. Since the observation time is 
much longer than the local-equilibrium time (125), it is appropriate to introduce the average 
free energy 

\ ptvar 



tvar 



F{t) dt , (127) 



averaging the local-equilibrium free energy (126) over the characteristic variation time of the 
alternating potential V{r,t). 

Suppose the alternating potential has the form 

V{r,t) = V{r)f{t), (128) 
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in which the temporal function f[t) varies in the interval 

fmin — fif) — fmax ■ (-^-^9) 

Let us define the function t — t{^) by the equation 

mo) = i ■ (130) 

By introducing the effective Hamiltonian 

H^ffiO = ^(^(0) , (131) 

we get 

HeffiO = H + j m^V{v),p{v) dv . (132) 
Also, we define the effective temperature 

^ mo) ^ . (133) 

The quantity 



plays the role of the distribution of the variable ^. With these notations, the free energy 
(127) transforms into 

f fmax 

F = - / %ff{0 In Tr exp {-/3e//(0^e//(0} P(0 ■ (135) 

Jmin 

Let ^* be a value in the interval fmin < ^* < fmax, and let us define 

T* = nffiC) = ^ . (136) 

Then, by the theorem of mean, the free energy (135) can be approximately represented as 

rfn 



F = —T' 



/Jmax 
In IV exp{-f3*Heff{0} p(0 ■ (137) 
J min 



This expression looks analogously to the free energy of a system in an external random 
potential ^V{r) that enters Eq. (132). 

To exemplify explicitly the form of distribution (134), let us take the time-dependent 
factor of the alternating potential (128) as 

f{t) = cos{u;t) . (138) 

The characteristic variation time here is, clearly, the period tyar — 27r/a;. The variation 
range of function (138), as defined in Eq. (129), is given by fmin — ~1 and fmax — 1- The 
function t{^) is found from Eq. (130), yielding 

t(^) = - arccosC . (139) 
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^Prom here, for distribution (134), we get 



1 



(-l<e<l)- 



(140) 



The variable ^, with distribution (140), is zero-centered, since 



and its dispersion is 



In this way, when considering the time-averaged behavior of a Bose system, subject to 
the action of an external alternating potential, we come to the picture that is similar to the 
description of a Bose system in an external spatially random potential. The consequences, 
therefore, should also be similar. The overall behavior of trapped atoms, in the presence of 
an alternating potential V{r,t), should be as follows. 

Let Vo be the strength of the alternating potential V{r.t) ~ Vq. The latter plays the 
role of Vji in the localization length (84). If the whole trap is perturbed by the alternating 
potential, then the characteristic trap length Iq = l/^/muJo, where is the trap frequency, 
plays the role of Ir in Eq. (84). Hence, for the localization length (84), we have 



When the amplitude of the alternating potential is small, such that Vq <^ ujq, then there is an 
extended condensate filling the trap. If the alternating potential oscillates with a frequency 
far detuned from any transition frequency for trapped atoms, then the extended condensate 
is just a perturbed ground-state condensate. But if the frequency of the oscillating potential 
is in resonance with the transition frequency between the ground-state condensate and an 
excited topological mode, then the extended condensate is formed by the fluctuating ground- 
state condensate mode and an excited coherent mode [61-64]. 

Increasing the amplitude Vq of the alternating potential perturbs the Bose-condensed 
system stronger. When Vo reaches the trap frequency loq, the locahzation length becomes 
comparable with the effective trap length Iq. At this moment, the granular condensate 
starts being formed. Since the applied external field is time dependent, the created granular 
condensate is nonequilibrium. This means that the disconnected regions with the condensate, 
with time, change their shapes, disappear and appear again, so that on average there always 
exist several such regions of the granular condensate. 

The nonequilibrium granular condensate can exist for the amplitude of the alternating 
potential in the range 



If the amplitude is so large that Vq ~ Uoylo/a, then lioc ~ a, and the granular condensate is 
completely destroyed. Then the whole system is a strongly nonequilibrium normal matter, 
containing no condensate, but being rather in a normal turbulent state. 
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The phase portrait of the nonequihbrium system, at zero temperature, quahtatively 
corresponds to that of the equihbrium random system, as is shown in Fig. 1, with the 
following analogies between the systems: 

equilibrium extended condensate nonequilihrium extended condensate, 

equilibrium granular condensate nonequilihrium granular condensate, 

equilibrium normal glass turbulent normal fluid. 

The principal difference between the equilibrium and nonequilibrium granular conden- 
sates is as follows. The equilibrium granular condensate represents a system with randomly 
distributed in space regions of the condensate, surrounded by the normal phase with no 
condensate, these regions being stationary in time, their location in space being fixed. The 
nonequilibrium granular condensate, at each instant of time, is analogous to its equilibrium 
counterpart. However, the regions with the condensate are not fixed in time. But, as time 
varies, the condensate regions change their shapes and locations. They can appear in new 
spatial locations, but disappear in others. What is fixed for a nonequilibrium granular con- 
densate is the average (over space and time) concentration of condensed atoms, so that the 
average condensate fraction Uq remains a well defined order parameter. At each instant 
of time, the system with the nonequilibrium granular condensate is a kind of a snapshot, 
which is analogous to the equilibrium granular condensate. Therefore, if trapped atoms are 
released from the trap, the time-of-flight observations for a nonequilibrium system will be 
similar to those for an equilibrium system. 

We may note that the nature of bosons, considered above, can be any. These could be 
usual bosonic atoms [2-5]. Or these could be composite bosonic molecules formed of fermions 
[3,65]. 

The discussed analogies between the random Bose-condensed systems and the nonequlib- 
rium Bose systems can be tested experimentally. Such experiments are now in progress in 
the Institute of Physics of Sao Carlos, University of Sao Paulo, Brazil. Preliminary results 
confirm these analogies. But a detailed exposition of experiments will be done in separate 
publications. 
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Figure 1: Qualitative phase portrait of a Bose system in a spatially random potential at zero 
temperature on the plane of the gas parameter 7 and disorder parameter (. 
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